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摘要 : 本 文 利 用 激 波 管 测量 了 乙 烷 在 不 同 压 力 、 温 度 和 当量 比 下 的 着 火 延迟 期 ， 并 利用 CHEMKIN 软件 结 
合 GRI-Mech 3.0 机 理 对 预 混合 气 着 火 过 程 进行 了 模拟 计算 和 化 学 反应 动力 学 分 析 。 结 果 表 明 ,GRI-Mech 3.0 
机 理 不 能 很 好 的 预测 乙 烷 的 着 火 延 迟 期 。 通 过 对 着 火 过 程 进行 敏感 性 分 析 ， 得 出 反应 CHs+02=C2Hs+HO， 
和 C2H4+H+M=C2H5s+M 对 着 火 延迟 期 的 影响 较 大 ， 并 对 其 反应 速率 常数 进行 了 修正 。 修 正 后 的 机 理 能 够 较 
好 地 改善 对 乙 烷 着 火 延迟 期 的 预测 ， 并 同时 保持 了 对 层 流 燃烧 速度 的 准确 预测 。 
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Study on the chemical kinetic mechanisms of ethane ignition 


CHEN Yi-Zhen, HU Er-J iang， ZHANG Zi-Hang, PAN Lun, HUANG Zuo-Hua 
(State Key Laboratory of Multiphase Flow in Power Engineering, Xi'an Jiaotong University, Xi’an 710049, 
China) 
Abstract: Ignition delay times of ethane were measured behind reflected shock waves under 
different pressures, temperatures, and equivalence ratios. Simulations and chemical kinetic 
analysis on the ignition process of the mixtures were carried out using CHEMKIN II package with 
GRI-Mech 3.0 mechanism.The simulation results showed large scatters with the experimental data 
especially at lean mixtures. Reaction C2?H5s+O2=C2H4+HO， and reaction 
CHst+H+M=CHs+Mwere found to havesignificantinfluenceon the ignition process of ethane, and 
their rate constants were modified in order to give more accurate prediction. The optimized 
mechanism can give better prediction of ignition delay times of ethane and maintains its good 
prediction onthe laminar flame speed. 
Key words:ethane; ignition delay time; GRI-Mech 3.0; chemical kinetic; sensitivity analysis 
0 下 过 学 反应 机 理 ， 都 是 十 分 重要 的 。 

着 火 延 迟 期 是 描述 燃料 着 火 特性 的 基础 参数 ， 
当今 世界 面临 着 化 石 能 源 不 可 再 生 并 日 益 枯 竭 。 是 实际 燃烧 装置 燃烧 室 设 计 和 优化 的 重要 依据 ， 同 
的 问题 ， 同 时 污染 气体 排放 与 温室 气体 控制 的 压力 ” 时 也 是 构建 燃料 动力 学 模型 的 必要 数据 。 国 际 上 对 
巨大 ， 因 此 提高 化 石 能 源 的 燃烧 效率 并 同时 寻找 清 于 乙 烷 着 火 延 迟 期 的 测量 和 研究 已 经 有 很 多 。 在 
洁 的 蔡 代 燃料 是 全 球 共 同 的 研究 热点 。 有 关 低 碳 燃 1970 年 , Bowman] 率 先 使 用 激 波 管 测量 了 乙 烷 在 压 
料 的 燃烧 基础 研究 在 解决 这 两 大 难题 的 过 程 中 占有 ”” 力 为 0.2-0.8 MPa 下 的 着 火 延 迟 期 。Cooke 等 人 让 在 
十 分 重要 的 地 位 ， 其 中 就 包括 乙 烷 的 燃烧 与 化 学 反 ”1975 年 测量 了 乙 烷 在 高 温 下 1400-2200 K 的 着 火 延 
应 动力 学 研究 。 首 先 ， 乙 烷 是 天 然 气 中 次 于 甲烷 的 。 述 期 。Hidaka 等 人 P11 利用 激 波 管 研究 了 乙 烷 的 热 解 
最 重要 的 一 种 组 分 ， 而 天 然 气 则 是 目前 发 展 最 为 成 和 氧化 机 理 并 测量 了 其 中 间 物 种 的 浓度 。Shim 等 人 
功 的 代用 燃料 之 一 。 其 次 ， 乙 烷 作为 中 间 产 物 广 泛 。 ”中 和 Lamoureux 等 人 四 也 都 在 激 波 管 上 测量 过 乙 烷 的 
地 存在 于 碳 氧化 物 的 氧化 过 程 中 ， 而 且 许多 较 大 分 。 着 火 延 迟 期 -Tranter 等 人 中 利用 激 波 管 研究 了 超 高 压 
子 碳 氧 化 物 的 化 学 反应 动力 学 机 理 都 是 以 低 碳 碳 氧 力 下 (34 MPaand 61.3 MPa) 乙 烷 的 热 解 和 氧化 ， 发 
化 物 的 化 学 反应 动力 学 机 理 为 基础 构建 的 。 因 此 ， 现 已 有 机 理 能 较 好 预测 超 高 压力 下 乙 烷 着 火 延 迟 期 。 
详细 而 准确 地 研究 乙 烷 的 化 学 反应 动力 学 ， 对 于 理 de Vries 等 人 由 测量 了 乙 烷 的 着 火 延迟 期 并 与 前 人 的 
解 乙 烷 的 燃烧 氧化 过 程 ， 和 发 展 更 大 分 子 燃料 的 化 ”实验 结果 作 了 对 比 ， 发 现 不 同 研究 者 测 得 的 乙 烷 着 
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火 延迟 期 的 拟 合 式 有 一 定 差 异 。 由 此 可 见 ， 虽 然 乙 
烷 的 着 火 延 迟 期 已 被 广泛 地 测量 过 ， 其 热 解 和 氧化 
机 理 也 已 经 被 较 详 尽 地 研究 ， 但 是 实验 数据 主要 集 
中 于 压力 在 0.02-1.0MPa， 当 量 比 为 0.5-2.0 之 间 。 在 
更 广泛 工 况 下 更 丰富 的 实验 数据 有 利于 乙 烷 着 火化 
学 反应 动力 学 机 理 进一步 的 发 展 。 

本 文采 用 激 波 管 和 CHEMKIN 软件 研究 乙 烷 在 
不 同 温度 和 压力 下 的 着 火 延 迟 期 。 研 究 混合 气 初始 
参数 〈 温 度 、 初 始 压 力 和 燃料 浓度 ) 对 着 火 延迟 期 
的 影响 ， 重 点 研究 和 验证 了 乙 烷 的 化 学 反应 动力 学 
机 理 一 一 GRI-Mech 3.0 机 理 。 该 研究 获得 的 基础 实 
验 数 据 既 可 以 丰富 乙 烷 基础 燃烧 数据 ， 同 时 也 对 发 
展 乙 烷 化 学 动力 学 模型 起 到 促进 作用 。 


1 实验 方法 与 计算 方法 


1.1 实验 : 
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图 1 着 火 延 迟 期 定义 


Fig. 1 The definition of ignition delay time 
图 1 给 出 了 本 文 着 火 延 迟 期 的 定义 。 着 火 延迟 


期 定义 为 入 射 激 波 到 达 激 波 管 端面 的 时 刻 与 CH* 信 


装置 和 方法 
本 文 测量 着 火 延迟 期 的 实验 利用 激 波 管 实验 平 
台 完 成 。 由 不 锈 钢 钢 管 组 成 的 激 波 管 管 体 被 PVC 膜 
分 成 低 、 高 压 段 两 部 分 , 长 度 分 别 为 5.3 m 和 2.0 m。 
激 波 管 在 实验 前 抽 真 空 至 1x103kPa 以 下 。 三 个 压力 


号 曲线 最 大 仙 


隔 。 首 先 ， 通 过 


寸 位 于 璧 


| 率 处 切线 在 基线 处 交点 之 间 的 时 间 间 
面 及 端面 的 四 个 压力 传感器 


(PCB, B26) 测 得 压力 数据 和 压力 升 高 的 时 刻 , 并 触发 
三 台 位 于 壁面 的 时 间 间 隔 计 数 器 (ECLUKE PM6690)。 


传感器 (PCB113B26) 安装 在 距离 实验 段 端面 不 同 
立 置 处 ， 用 于 测量 入 射 激 波 速 度 ， 并 通过 一 维 理想 
ea dar tenn 实验 混 

气 在 不 锈 钢 混 气 灌 中 混合 24 h 以 上 ， 配 气 前 混 气 
维 抽 真空 至 1x103kpa 以 下 。 通过 安装 在 低压 
处 的 压力 传感器 (PCB113B03) 和 相同 位 置 的 光电 
倍增 管 (PMT ) 共同 诊断 混合 气 的 着 火 。PMT 前 放 


置 窄带 通 滤 光 片 , 用 于 探测 波长 430nm 的 CH* 自 发 
光 信 号 。 详 细 的 实验 装置 介绍 见 文献 四。 
表 1 激 波 管 实 验 混合 气 组 分 
Table 1 Compositions of the test mixtures 
混合 气 。 $$ Xpand% Toxy% XW p/MPa 
1 0.5 1.0 7.0 92.0 0.12, 0.5, 2 
2 1.0 1.0 3.5 95.5 0.12, 0.5, 2 
3 2.0 1.0 1.75 97.25 0.12, 0.5, 2 
4 1.0 0.75 2.625 96.625 0.5 
5 1.0 1.5 5.25 93.25 0.5 
6 1.0 2.0 7.0 91.0 0.2,0.5 
7 0.5 0.25 1.75 98.0 0.11 
实验 所 用 的 乙 烷 纯度 为 99.9%，0, 纯度 为 
99.995%，Ar 和 He 纯度 为 99.999%。 实 验 工 况 为 : 


当量 比 为 0.5、1.0、2.0。 着 火 压力 p=0.12-2.0 MPa， 
混合 气 中 各 成 分 的 摩尔 分 数 如 表 1 所 示 。 


然后 ， 对 所 得 的 三 个 时 间 间 隔 进行 线性 插值 可 


以 推 


Para 


管 测 得 。 


算 入 射 激 波 速度 .CH* 信 号 是 由 位 于 端 


的 光电 倍增 


1.2 计算 方法 和 化 学 反应 机 理 
本 文 还 利用 化 学 反应 动力 学 软件 CHEMKIN 中 


的 SENKIN 模块 对 乙 烷 的 


自 着 火 特 性 进行 了 数值 模 


拟 和 动力 学 


分 析 。 本 文 所 用 机 理 为 GRI-Mech 3.0 机 


性 分 析 并 进行 修正 ， 修 1 


理 bl. 同 时 ， 基 于 本 文 的 实验 数据 和 DeSain 等 人 0 
的 计算 结果 ， 本 文 对 GRI-Mech 3.0 机 理 进 行 了 敏感 


E 后 的 机 理 


较 好 地 


能 够 同时 


证 


预测 乙 烷 的 着 火 延 j 
结果 与 分 析 


迟 期 和 层 流 燃 烧 速 度 。 


2.1 着 火 延 迟 期 实验 值 的 验证 


图 2 给 出 了 在 压力 为 0.11 MPa， 当 量 比 为 0.5， 


乙 烷 摩尔 分 数 为 外 mawe=0.25% 时 ,以 及 压力 为 0.2MPa， 
当量 比 为 1， 乙 烷 摩 尔 分 数 为 系 rane=2% 时 ，deVries 


等 四 的 实验 结果 与 本 下 


在 这 两 个 工 况 下 本 文 的 实验 结果 与 deVries 等 
合 很 好 ， 说 明 本 文 的 实验 结果 
乙 烷 的 化 学 反应 动力 学 机 理 有 很 


验 结果 
重复 性 ， 对 于 验证 
好 的 参考 。 


究 的 实验 结果 对 比 。 可 以 看 出 ， 


人 的 实 
\ 有 较 好 的 
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2 着 火 延迟 期 实验 值 对 比 


Fig. 2 Comparison between present measured data and the 


measured data from reference 


2.2 着 火 延 迟 期 实验 值 与 计算 值 对 比 


图 3 给 出 了 乙 烷 着 火 延迟 期 实验 测量 值 和 计 和 介 
值 的 对 比 。 点 为 激 波 管 实验 值 , 虚线 为 GRI-Mech 3.0 
机 理 的 计算 值 ， 实 线 为 修正 后 GRI-Mech 3.0 机 理 的 
计算 值 。 由 图 可 知 ， 在 当量 比 和 燃料 浓度 不 变 的 情 
况 下 ， 初 始 压 力 越 高 乙 烷 的 着 火 延 迟 期 越 短 ， 而 在 
初始 压力 和 当量 比 不 变 的 情况 下 ， 人 燃料 浓度 越 高 乙 
烷 的 着 火 延迟 期 越 得 ， 说 明 初 始 压 力 和 燃料 浓度 的 
提高 对 乙 烷 的 着 火 有 促进 作用 。 由 乙 烷 着 火 延迟 期 
的 计算 值 可 以 看 出 ,在 当量 比 为 0.5 和 1.0 的 工 况 下 ， 
GRI-Mech 3.0 对 乙 烷 着 火 延迟 期 的 计算 值 都 明显 高 
于 实验 值 ， 在 个 别 工 况 下 两 者 的 差别 甚至 达到 一 个 
数量 级 。 因 此 ， 本 文 对 GRI-Mech 3.0 机 理 进行 了 修 
正 ， 在 这 两 种 工 况 下 ， 而 修正 后 的 机 理 能 够 较 好 二 
预测 乙 烷 的 着 火 延 迟 期 。 而 在 当量 比 为 2.0 的 工 况 1 
原 机 理 与 修正 后 机 理 的 预测 情况 均 不 理想 ， 有 待 进 

步 的 研究 。 以 下 部 分 通过 化 学 反应 动力 学 分 析 详 


细 地 图 述 了 对 机 理 的 修正 过 程 。 
2.3 化 学 反应 动力 学 分 析 


为 了 确定 各 基 元 反应 对 着 火 延 迟 期 的 影响 ， 本 
文 对 乙 烷 的 着 火 延迟 进行 敏感 性 分 析 。 敏 感性 分 析 
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图 3 着 火 延迟 期 实验 值 与 计算 值 对 比 〈 点 : 实验 值 ， 虚 


: GRI-Mech 3.0 计算 值 ， 实 线 : 修正 后 的 GRI-Mech 3.0 计 


算 值 ) 


Fig. 3Comparison between experimental and modeling 
results (Symbols: measurement; dash dot lines: GRI-Mech 3.0; 
solid lines: Modified GRI-Mech 3.0) 


是 指定 量 地 研究 某 一 因素 发 生变 化 时 对 结果 的 影响 
程度 ， 根 据 Petersen 等 人 0 的 定义 ， 着 火 延迟 期 的 
敏感 性 系数 如 下 公式 所 示 : 


r(2. 0k,) — zr(0. 5k,) 
1. 57(k,) 


其 中 7 为 着 火 延迟 期 ,为 第 i 个 基 元 反应 速率 系数 。 
敏感 性 系数 为 正 的 反应 生成 稳定 的 产物 或 自由 基 
降低 总 体 反 应 活性 ， 使 着 火 延 迟 期 增加 ， 而 敏感 性 
系数 为 负 的 反应 能 够 生成 活性 更 强 的 自由 基 或 更 多 
的 自由 基 , 增强 总 体 反应 活性 , 使 着 火 延迟 期 缩短 。 
由 于 低压 时 原 机 理 的 着 火 延 迟 预 测 偏差 较 大 ， 故 主 
要 对 该 工 况 进行 敏感 性 分 析 。 图 4 给 出 了 p=0.12MPa,， 
7=1250 K，Y=0.5 条 件 下 GRI-Mech 3.0 机 理 着 火 延 
述 期 的 敏感 性 系数 。 与 其 他 碳 氧 燃料 相似 ， 在 高 温 
条 件 下 R8 (HH+HO=O+OH) 是 敏感 度 最 大 的 促进 着 
火 的 反应 。 在 列 出 的 敏感 性 系数 最 大 的 11 个 反应 中 ， 
与 乙 烷 燃 料 分 子 相 关 的 反应 有 6 个 : 

R74: CH4+HCHMD)=C2Hs(CHMD 

R78: CHer+H=C2Hs+H2， 

R112: C2H4+OH=C2H3+H2O 

R158: 2CH3(+MD)=C2He(HMD) 

R175: C"2Hs+O2=C2H4+HO。， 

R294: C2Hi+O2=CHCHO+O 


其 余 均 为 小 分 子 物种 反应 ， 这 些小 分 子 物种 反应 在 
验证 其 它 燃 料 如 氢气、 甲烷 等 燃烧 特性 时 ， 已 经 被 
大 量 的 实验 结果 所 验证 ， 故 本 文 对 GRI-Mech 3.0 机 
理 的 修正 主要 集中 在 这 6 个 与 乙 烷 燃料 分 子 相 关 的 
反应 上 。 


Se 


(1) 


R38: H+O2=O+OH 
R74: H+C2H4(+M)=C2HS(+M) 
R158: 2CH3(+M)=C2HE(+M) 
R119: HO2+CH3=OH+CH3O 
R116: 2HO2=O2+H2O2 

R112: OH+C2H4=C2H3+H2O 
R294: C2H3+02=O+CH2CHO 
R87: OH+HO2=O2+H2O 

R45: H+HO2=O2+H2 

R78: C2H6+H=C2H5+H2 


R175:; C2H5+O2=HO2+C2H4 


图 4 着 火 延迟 期 的 敏感 性 系数 
Fig.4Sensitivity analysis of ignition delays 
Zhang 等 人 号] 在 对 乙 烷 进行 反应 路 径 分 析 时 发 
现 乙 烷 大 部 分 是 被 生成 乙 基 的 脱氧 反应 所 消耗 。 生 
成 的 乙 基 会 进一步 发 生 脱 氨 反 应， 其 中 部 分 乙 基 与 
氧 原 子 发 生 反 应 生成 乙烯 和 HO; 自由 基 。DeSain 等 


人 0 在 2003 年 提出 了 C2Hs+O; 的 反应 速率 常数 的 理 
论 计算 值 ， 它 包含 多 个 反应 通道 。 由 于 本 文 的 研究 
范围 主要 集中 于 高 温 区 , 而 在 高 温 区 C2H5+0O2 的 反 
应 通道 主要 是 CsH; 与 0, 的 直接 脱毛 反应 ， 即 
CHs+Os=C2Hst+HO,, 因此 在 GRI-Mech 3.0 机 理 中 只 
考虑 这 一 反应 通道 具备 一 定 的 合理 性 。 但 是 
GRI-Mech 3.0 中 R175 的 反应 速率 常数 的 值 与 DeSain 
等 人 的 理论 计算 值 相 差 很 大 。 在 将 CsH; 与 0; 的 直接 
脱氧 反应 的 DeSain 等 人 的 理论 计算 值 代 入 
GRI-Mech 3.0 后 ， 预 测 的 着 火 延 迟 期 明显 降低 ， 但 
是 在 个 别 工 况 下 仍然 高 于 实验 值 。 此 外 ， 反 应 中 大 
部 分 乙 基 则 会 发 生 单 分 子 反应 生成 乙烯 与 H 自由 基 ， 
从 而 增加 反应 中 H 自由 基 的 浓度 。 事 实 上 链 分 支 反 
应 R74 能 明显 促进 乙 烷 的 着 火 。 考 虑 到 原 机 理 预 测 
的 乙 烷 着 火 延 迟 期 过 长 ， 在 将 R74 的 反应 速率 常数 
的 高 压 极限 和 低压 极限 同时 提高 一 倍 以 后 ， 修 正 
的 机 理 能 够 相对 较 好 地 预测 着 火 延 迟 期 的 变化 。 本 
文 的 实验 与 模拟 工作 为 R74 反应 速率 常数 的 改进 提 
供 参 考 与 依据 ， 为 R74 的 反应 速率 常数 的 高 精度 计 
算 提 供 了 重要 参考 。 
为 了 进一步 验证 修正 后 机 理 对 乙 烷 层 流 燃 烧 速 
度 的 预测 ， 图 5 给 出 了 GRI-Mech 3.0 对 乙 烷 的 层 流 
燃烧 速度 的 计算 值 与 修正 后 机 理 的 计算 值 的 对 比 。 
从 图 中 可 以 看 出 ， 修 正 后 的 机 理 与 原 机 理 在 乙 烷 层 


@ Hassan MI. et al. [14] 
0 人 Jomaas G. et al. [15] 


Flame speed/ cmis 
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图 5 层 流 燃 烧 速 度 实验 值 与 计算 值 的 对 比 (点 : 实验 值 ， 虚 
线 : GRI-Mech 3.0 计算 值 ， 实 线 : 修正 后 的 GRI-Mech 3.0 计 
算 值 ) 


Fig. $ Comparison between experimental and modeling results of 


laminar flame speed (Symbols: measurement; dash dot lines: 
GRI-Mech 3.0; solid lines: Modified GRI-Mech 3.0) 
流 燃 烧 速 度 的 预测 上 差别 并 不 大 。 根 据 Hassan 等 人 
0 及 Jomaas 等 人 0 的 实验 值 , GRI-Mech 3.0 对 乙 烷 


的 层 流 燃 烧 速 度 的 预测 较为 准确 。 所 以 修正 后 的 机 
理 在 维持 了 原 机 理 对 乙 烷 层 流 燃 烧 速 度 的 较 准确 的 
预测 结果 的 基础 上 ， 对 GRI-Mech 3.0 的 乙 烷 着 火 延 
迟 期 的 预测 有 所 改善 。 


3 结 i Tb 


本 文 利用 激 波 管 和 化 学 反应 动力 学 软件 
CHEMKIN 开展 了 不 同 温度 、 压力、 当量 比 和 燃料 浓 
度 时 乙 烷 的 着 火 延 迟 期 测量 和 动力 学 分 析 ， 获 得 的 

结论 如 下 : 

(1) 获得 了 乙 烷 不 同 条 件 下 的 着 火 延 迟 期 ， 在 本 研 
究 工 况 范围 内 ， 压 力 的 增加 和 当量 比 的 增加 对 着 火 
有 明显 的 促进 作用 。 着 火 延迟 期 的 对 数 与 温度 的 倒 
数 呈 良好 的 线性 关系 。 
(2) 着 火 延 迟 期 的 实验 值 和 计算 值 对 比 发 现 ， 
GRI-Mech 3.0 机 理 预测 的 着 火 延 迟 期 偏 高 。 

《3 ) 通过 敏感 性 分 析 , 找 出 了 GRI-Mech 3.0 机 理 : 
影响 乙 烷 着 火 延 迟 期 最 主要 基 元 反应 ， 并 对 其 速率 


dT 


it 


TT 


常数 进行 了 优化 ， 修 正 后 的 机 理 能 够 较 好 的 预测 实 
验 结果 。 
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